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We present a first-principles investigation of structural, electronic and optical properties of bulk 
crystalline AgsN, AgN and AgN2 based on density functional theory (DFT) and many-body pertur- 
bation theory. The equation of state (EOS), energy-optimized geometries, cohesive and formation 
energies, and bulk modulus and its pressure derivative of these three stoichiometrics in a set of 
twenty different structures have been studied. Band diagrams and total and orbital-resolved density 
of states (DOS) of the most stable phases have been carefully examined. Within the random-phase 
approximation (RPA) to the dielectric tensor, the single-particle spectra of the quasi electrons and 
quasi holes were obtained via the GW approximation to the self-energy operator, and optical spec- 
tra were calculated. The results obtained were compared with experiment and with previously 
performed calculations. 
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I. INTRODUCTION 

It is well-known now that late transition-metal nitrides 
(TMNs) usually possess interesting properties leading 
to a variety of potential technological applications^"'^. 
Hence, a significant number of quantum mechanical ab 
initio calculations of the structural and physical prop- 
erties of this family of materials have appeared in the 
literature. 

Since Juza and Hahn^ succeeded to synthesize CU3N in 
1939, copper nitrides have been produced through vari- 
ous techniques and their properties and applications have 
been the subject of many theoretical and experimental 



published works^. Due to its early discovery, copper ni- 
tride may now be considered as the most investigated 
among the late TMNs*^. 

On the other hand, the nitride of silver, the next el- 
ement to copper in group 11 of the periodic table, has 
been known for more than two centuries^'*. However, de- 
spite its earlier discovery, silver nitride may be the least 
theoretically studied solid in the late TMNs family. Ex- 
perimental efforts to investigate structural*'^, electronic* 
and formation^'^"^^ properties of silver nitrides have been 
made by some researchers. 

In 1949, Hahn and Gilbert^ carried out the first* struc- 
tural study on the reported stoichiometry, AgaN. They 
claimed an fee structure with a ~ 4.369 A and Z ~ 4/3 
(i.e. 4 Ag atoms in the unit cell). A long time later in 
1982, Haisa^^ suggested that the Ag atoms are located 
at the corners and face centers of the unit cell, while the 
N atoms, which may be statistically distributed in the 
octahedral interstices, were given no definite positions^^. 

According to the calculated N radius, AgsN can be 
described as an ionic compound, and recent ab initio 
calculations on the proposed structure revealed insulat- 
ing characteristics with a fundamental band gap close to 
1.35 eV. On the other hand, due to the similar lattice of 
the parent Ag and the easily separated N as N2, it can 
also be argued that this compound is a metal, supporting 
its black color*. 

Under ordinary conditions^, it was found that sil- 
ver can form AgaN^'^ from ammoniacal solutions of sil- 
ver oxide *"'*. The black metallic- looking solid outcome, 
AgaN, is an extremely sensitive explosive compound 
It may explode due to the slightest touch, even from the 
impact of a falling water droplet but it is relatively 
easy to handle under water or ethanol*. The explosive 
power is due to the energy released during the decompo- 
sition reaction: 



2Ag3N — > 6Ag + N2 



(1) 



Even in storage at room temperature, this solid com- 
pound decomposes slowly according to Eq. 1 above*' 
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From a thermochemical point of view, it was found that 
there is no stable intermediate stage in this decomposi- 
tion, but there may be a metastable intermediate species 
(phase) with a remarkably low decomposition rate*". At 
this point, it may be worth mentioning that the ther- 
mochemistry of silver nitride systems is not fully docu- 
mented in standard handbook data*". 

In their 1991 work. Shanlcy and Ennis^ stated: "Many 
of the samples ... did not survive the minimum han- 
dling required to move them, container and all, to the 
X-ray stage. ... More vigorously explosive samples prop- 
agated throughout their mass leaving no visible residue. 
Even among supersensitive materials, silver nitride is a 
striking example of a compound "teetering on the edge 
of existence" . Under the circumstances, we did not suc- 
ceed in developing data on the proportion of silver nitride 
required for explosive behavior in these mixtures. " 

Thus, beside the potential hazard to lab workers due to 
its sensitive explosive behavior, characterization of silver 
nitride is hindered by its extremely unstable (endother- 
mic) nature , and we are presented with an incomplete 
picture of structural, electronic and optical properties of 
this material. Surprisingly, this lack of detailed knowl- 
edge of many physical properties of silver nitride stimu- 
lated only very few published ab initio studies. 

In the present work, first-principles calculations were 
carried out to investigate the lattice parameters, equation 
of state, relative stabilities, phase transitions, electronic 
and optical properties of silver nitrides in three different 
chemical formulae and in various crystal structures. Cal- 
culation methods are described in Sec. II. In Sec. Ill, 
results are presented, discussed and compared with ex- 
periment and with previous calculations. The article is 
concluded with some remarks in Sec. IV. 



II. CALCULATION METHODS 
A. Stoichiometries and Crystal Structures 

To the best of our knowledge, the only experimen- 
tally reported stoichiometrics of Ag-N compounds are 
AgsN'' and AgNs*". However, previous ab initio stud- 
ies on Ag-N compounds considered Ag4N^, AgsN^'^^, 
Ag2Ni^, AgNi^'^'^ and AgN2^^'^'' in some cubic struc- 
tures only. Consideration of stoichiometries other than 
the reported ones is probably due to the fact that many 
transition metals nitrides (TMs) are known to form more 
than one nitride^*. Hence, our interest in investigating 
AgN and AgN2 is based on this fact. 

For AgsN, we consider the following seven structures: 
the face-centered cubic structure of AlFca (DO3), the sim- 
ple cubic structure of CraSi (A15), the simple cubic struc- 
ture of the anti-ReOs (DO9), the simple cubic structure 
of AgaAu (LI2), the body-centered cubic structure of 
C0AS3 (DO2), the hexagonal structure of e-FesN, and the 
trigonal (rhombohedric) structure of RhFs. 

For AgN, the following four structures were considered: 



the face-centered cubic structure of NaCl (Bl), the sim- 
ple cubic structure of CsCl (B2), the face-centered cubic 
structure of ZnS zincblende (B3), the hexagonal structure 
of NiAs (B81), the hexagonal structure of BN (Bk), the 
hexagonal structure of WC (Bh), the hexagonal structure 
of ZnS wurtzite (B4), the simple tetragonal structure of 
PtS coopcrite (B17), and the face-centered orthorhombic 
structure of TIF (B24). 

AgN2 was studied in the following nine structures: the 
face-centered cubic structure of CaF2 fluorite (CI), the 
simple cubic structure of FeS2 pyrite (C2), the simple 
orthorhombic structure of FeS2 marcasite (C18) and the 
simple monocline structure of CoSb2 (CoSb2). 



B. Electronic Relaxation Details 

In this work, electronic structure spin density func- 
tional theory (SDFT)^^'^° calculations as implemented 
in the VASP code^^"^^ have been employed. To self- 
consistently solve the Kohn-Sham (KS) equations^^ 

j-^V^+Z^r'^-l-Wr) 
2me J |r - r'l 

^ . (2) 

+F^;^[n(r)]Ur^r)=,r''V'r''(r), 



where i, k and a are the band, k-point and spin in- 
dices, receptively, VASP expands the pseudo part of the 
KS one-particle spin orbitals ip^'^lr) on a basis set of 
plane- waves (PWs). Only those PWs with cut-off en- 
ergy Ecut < 600 eV have been included. The Brillouin 
zones were sampled using F-centered Monkhorst-Pack^^ 
17x17x17 meshes. Any increase in the Ecut value or 
in the density of the k-mesh produces a change in the 
total energy less than 3 mel^/atom. For static calcula- 
tions, partial occupancies were set using the tetrahedron 
method with Blochl corrections^^"'^^, while the smearing 
method of Methfessel-Paxton (MP)'^^ was used in the 
ionic relaxation, and Fermi surface of the metallic phases 
has been carefully treated. The Perdew-Burke-Ernzerhof 
(PBE) paramctrization'^'^"^^ of the generalized gradi- 
ent approximation (GGA)'^^"'^^ was employed for the 

exchange-correlation potentials yycl^l'^)]- T^^*^ imple- 
mented projector augmented wave (PAW) method^®''^^ 
was used to describe the core- valence interactions Vext{^), 
where the Ad^^Bs^ electrons of Ag and the 25^2^"^ elec- 
trons of N are treated explicitly as valence electrons. 
While for these valence electrons only scalar kinematic 
relativistic effects are incorporated, the PAW potential 
treats the core electrons in a fully relativistic fashion^^. 
No spin-orbit interaction of the valence electrons has 
been considered. 
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C. Geometry Relaxation and EOS 



E. GWA Calculations and Optical Properties 



At a set of externally imposed lattice constants, ions 
with free internal parameters were relaxed until all force 
components on each ion were less than 0.01 eV/ K. This is 
done following the implemented conjugate-gradient (CG) 
algorithm. After each ion relaxation step, static total en- 
ergy calculation with the tetrahedron method was per- 
formed, and the cohesive energy per atom Ecoh was cal- 
culated from 



Although a qualitative agreement between DFT- 
calculated optical properties and experiment is possible, 
accurate quantitative description requires treatments be- 
yond the level of DFT^^. Another approach provided 
by many-body perturbation theory (MBPT) leads to a 
system of quasi-particle (QP) equations, which can be 
written for a periodic crystal as^^"**^ 



E 



Ag„N„ 



E. 



Ag,„N,. 
solid 



Z X mE. 



Z X 



n) 



, (3) 



where Z is the number of AgmN„ formulas per unit cell, 
-^i^om *5.nd i?Som ^-re the energies of the isolated non- 
spherical spin-polarized atoms, and m,ri = 1,2 or 3 arc 
the stoichiometric weights. The obtained E^oh as a func- 
tion of volume V per atom were then fitted to a Birch- 
Murnaghan 3rd-order equation of state (EOS)^° and the 
equilibrium volume Vq, the equilibrium cohesive energy 
Eq, the equilibrium bulk modulus 



Bn = -V 



,dP 
dV 



V=Vo 



and its pressure derivative 



(4) 



(8) 



Practically, the wave functions tpf^{^) are taken from 
the DFT calculations. However, in consideration of 
computational cost, we used a less dense mesh of It- 
points (10 X 10 X 10). The term E(r, r'; e^^) is the self- 
energy which contains all the exchange and correlation 
effects, static and dynamic, including those neglected in 
our DFT-GGA reference system. In the so-called GW 
approximation^^, E is given in terms of the Green's func- 
tion G as 



° " dP p=o 



were determined. 
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V=Vo 



(5) 



D. Formation Energy 

Beside the cohesive energy, another measure of relative 
stability is the formation energy Ef. Assuming that sil- 
ver nitrides AgmN„ result from the interaction between 
the N2 gas and the solid Ag(Al) via the reaction (com- 
pare with Eq. 1) 



mAg' 
Ef can be given by 
EfiAg^K°''') = 



lid 



-Nr 



As N''' 



ilid 



(6) 



£;eoh(Ag„Nf''i) 
m£;coh(Ag^°'''^) + f£;coh(Nr) 



(7) 



Here m,n ~ 1,2,3 are the stoichiometric weights and 
£'coh(Ag,„N^°'"^) is the cohesive energy per atom as in Eq. 
3. The cohesive energy i?coh(Ag''°'"^) and other equilib- 
rium properties of the elemental metallic silver are given 
in Table I. The cohesive energy of the diatomic nitro- 
gen (i?coh(N|^'')) was found to be —5.196 eF/atom corre- 
sponding to an equilibrium N-N bond length of 1.113 A 
(For more details, see Ref. 5). 



I de'G(r,r';e,e')W^(r,r';e), 



(9) 



where the dynamically (frequency dependent) screened 
Coulomb interaction W is related to the bare Coulomb 
interaction v through 



W{r,r';e) = j J drie \r, n; e)?;(ri, r'). 



(10) 



with £, the dielectric matrix, calculated within the so- 
called random phase approximation (RPA). Wc followed 
the GWq self-consistent routine on G, in which the QP 
eigenvalues 



Re 



((*? 



QP 

k 



EtKS — Vxc + SgWo 



(11) 



are updated in the calculations of G, while W is kept 
at the DFT-RPA level. Four updates were performed, 
and after the final iteration in G, e is recalculated within 
the RPA using the updated QP eigenvalues''^'''*'*^. From 
the real £-[c{lij) and the imaginary £ini(a;) parts of this 
frequency-dependent microscopic dielectric tensor, one 
can derive all the other frequency-dependent dielectric 
response functions, such as reflectivity i? (w), transmitiv- 
ity T (w) = 1 — R{uj), refractive index n{uj), extinction 
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coefficient k (oj), and absorption coefficient a (w)^'' 

2 



[ere (t^) + jeim (w)]^ + 1 

1 / - 

« (^) = ( [£?c (^) + eL (W)] ^ + Ere (w) 



K (lu) 

a (lu) 



(13) 
(14) 
(15) 
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III. RESULTS AND DISCUSSION 
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The energy- volume equation of state (EOS) for tire dif- 
ferent structures of AgsN, AgN2 and AgN are depicted 
in Figs. 1, 2 and 3, respectively. The corresponding cal- 
culated equilibrium properties are given in Tabic I. In 
this table, we ordered the studied phases according to 
the increase in the nitrogen content; then within each se- 
ries, structures arc ordered in the direction of decreasing 
structural symmetry. For the sake of comparison, we also 
presented results from experiment and from previous ab 
initio calculations; and, whenever appropriate, the cal- 
culation method and the XC functional are also given in 
footnotes of the Table. 

The calculated equilibrium properties: cohesive 
energies, formation energies, volume per atom, volume 
per Ag atom, and bulk modulus and its pressure 
derivative which arc given Table I, arc visualized in 
Fig. 4. This kind of visualization allows us to study 
the effect of nitridation on the parent Ag(Al), since 
all quantities in this figure are given relative to the 
corresponding ones of the elemental Ag(Al) given in 
the first row of Table I. Moreover, one can easly com- 
pare the properties of these phases relative to each other. 



FIG. 1. (Color online.) Cohesive energy _Ecoh(eI^/atom) ver- 
sus atomic volume V (A"^/atom) for AgsN in seven diflferent 
structural phases. 




FIG. 2. (Color online.) Cohesive energy i5coh(ey/atom) ver- 
sus atomic volume V (A^/atom) for AgN in nine different 
structural phases. 



A. EOS and Relative Stabilities 

Considering -Ecoh in the AgaN series, Fig. 1 shows 
clearly that the E{V) relations of AgaN in DOg, DO2 and 
RhFa phases are almost identical, corresponding to equi- 
librium cohesive energy (Table I) of —2.513, —2.514 and 
—2.514 eV/ atom, respectively. This behavior in the EOS 
could be traced back to the structural relationships be- 
tween these three structures, since both DO2 and RhFa 
can simply be derived from the more symmetric DO9 (see 
Rcf. 5 and references therein). These structural relations 
may reflect in the EOS's and in other physical properties, 
and the three phases may co-exist during the AgsN syn- 
thesis process. Relative to the elemental Ag, these three 
phases tend not to change the -Bcoh (Fig- 4), lowering 
it only by 0.03 el^/atom, as can been seen from Ta- 
ble I. It may be worth to mention here that the simple 




FIG. 3. (Color online.) Cohesive energy E coh{eV / aiom) ver- 
sus atomic volume 1/ (A'' /atom) for AgN2 in four different 
structural phases. 
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TABLE I. Calculated and experimental zero-pressure properties of the twenty studied phases of AgaN, AgN and AgN2: Lattice 
constants (a(A), 6(A), c(A), a(°) and /3(°)), equilibrium atomic volume Vb(A''/atom), cohesive energy iJcoh(eV^/atom), bulk 
modulus Bo{GPa) and its pressure derivative Bq, and formation energy E f {eV / atom) . The presented data are of the current 
work (Pres.), experimentally reported (Expt.) and of previous calculations (Comp.). 



Structure 




a(A) 




6(A) 


c(A) 


a 


°) °' (3(°) 




(A'^ /atom 


) -Bj.„j^(eV7atoiii) 


B„(GPa) 




E 


f (e V'/ atom) 




Ag 




Pres. 


4.164 












18 


06 


— 2.484 


88.188 


5.793 






Al 


Expt. 


(4.0S57G ± 0.00018)^ 














-2.95'' 


100.7'=, 101° 


6.12'^ 








Comp. 


4.01* 
















-3.59S, -2.66'', 
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142' 


5.00^ 5.70j, 
5.97^ 






AggN 


D03 


Pres. 


6.322 












15 


79 


-2.055 


98.356 


5.457 


1 


107 


A15 


Pres. 


5.065 












16 


24 


-1.976 


92.280 


5.470 


1 


186 


DOg 


Pres. 
Comp. 


4.328 
3.995I 


4.169', 4.292^ 










20 


27 


-2.513 


71.980 
95.7', 87.1= 


5.386 





649 


L12 


Pres. 


3.972 












15 


67 


-2.081 


100.743 


5.530 


1 


081 


D02 


Pres. 


8.662 












20 


31 


— 2.514 


72.230 
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Pres. 


5.967 






5.560 






21 


43 


— 2.469 


64.737 


2.335 
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RhF3 


Pres. 
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= 59.989 


20 


31 


— 2.514 


72.237 


5.396 
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Expt. 


4.369" 
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AgN 
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4.617 
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-2.253 
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5.145 


1 
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Comp. 


4. 571, 4.506', 4.619*^ 
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4.476° 


4.606P 
















197.18°, 147. 40P 


4.883° 








Pres. 


2.873 












11 


86 


-2.021 


146.157 


5.260 


1 


819 


B2 


Comp. 


2.833"', 2.806', 2.876^ 
2.78°, 2.87P 
















138. 96P 
204.10° 


4. 8231^ 
5.451° 








Pres. 


4.950 












15 


16 


-2.122 


109.639 


5.210 


1 


718 


B3 


Comp. 


4.881, 4.816', 4.946'' 
4.79°, 4.94P 
















100. llP 
151.05° 


5.825f 
4.542° 






BSj 


Pres. 


3.544 






4.929 






13 


40 


-1.996 


130.485 


5.240 


1 


844 


Bu 


Pres. 


3.521 






9.368 






25 


15 


-1.891 


57.077 


5.110 


1 


949 


Bh 


Pres. 


3.096 






3.023 






12 


55 


-2.121 


141.385 


5.285 


1 


719 


B4 


Pres. 
Comp. 


3.501 
3.41°, 


3.54P 




5.734 

5.52°, 5.69 


P _ 




15 


22 


— 2.113 


105.992 

143.68°, 110. I2P 


5.467 

4.82°, 4.663P 


1 


727 


Bir 


Pres. 


Pres. 3.158 




5.560 






13 


86 


— 2.517 


132.556 


5.185 


1 


323 


B24 


Pres. 


4.337 




4.601 


5.091 






12 


70 


-2.202 


138.704 


5.132 


1 


638 


AgN2 




Pres. 


5.157 












11 


43 


-1.959 


164.844 


4.996 


2 


333 


CI 


Comp. 


5.124I, 5.055', 5.172^ 
5.013'", 5.141" 
















181.3', 164.5= 
215", 164" 








C2 


Pres. 


5.617 












14 


77 


-3.626 


30.058 


6.894 





666 


CIS 


Pres. 


3.440 




4.513 


5.508 






14 


25 


-3.680 


35.878 


7.269 





612 


CoSb2 


Pres. 


5.976 




5.651 


10.261 




= 151.225 


13 


90 


-3.699 


35.117 


7.822 
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vhilc bulk mudulii arc given at room temperature. 



^ Ref.^". Tills is an average of 56 experimental values, at 20° C. 
^ Ref.^^. Cohesive energies are given at X and 1 atm = 0.00010 GPa; 

Ref. (25) in 52: at room temperature, 
d See Refs. (8)-(ll) in 52. 

' Ref.^^. using the full-potential linearized augmented plane waves (LAPW) method within LDA. 
using projector augmented wave (PAW) method within LDA. 
using projector augmented wave (PAW) method within GGA(PW91). 
using projector augmented wave (PAW) method within GGA(PBE). 

using a scmiempirical estimate based on the calculation of the slope of the shock velocity vs. particle velocity curves obtained 
from the dynamic high-pressure experiments. The given values are estimated at ~ 298 K. 

Rcf.^^: using a scmiempirical method in which the experimental static P — V data are fitted to an EOS form where Bq and B'q are 
adjustable parameters. The given values are estimated at ~ 298 K. 

Rcf.^^: using the so-called method of transition metal pseudopotcntial theory; a modified form of a method proposed by Wills and 
Ifarrison to represent the effective interatomic interaction. 

using the full-potential linearized augmented plane waves (LAPW) method within LDA. 
using the full-potential linearized augmented plane waves (LAPW) method within GGA. 

using full-potential (linearized) augmented plane waves plus local orbitals (FP-LAPW+lo) method within LDA. 
using using full-potential (linearized) augmented plane waves plus local orbitals (FP-LAPW-l-lo) method within GGA(PBE). 
using pseudopotcntial (PP) method within LDA. 

using linear combinations of atomic orbitals (LCAO) method within LDA. Bq's are calculated from elastic constants, 
using linear combinations of atomic orbitals (LCAO) method within GGA. Bq's are calculated from elastic constants. 
This is the face centered cubic (fee) structure with Z = 4/3 (i.e. 4 Ag atoms in the unit cell) suggested by Hahn and Gilbert according 
to some measurements (Ref. 9). 
Ref 9. 
Ref 8. 

This is the average of the experimental values: (-1-314.4 T 2.5) kj/mol'^ = (3.25853 ± 0.02591) eV/atom, +199 kj/mol" = 
2.062 eV/atom, +255 kj/mol^ = 2.643 eV/atom, and +230 kj/mol^'' = 2.384 eF/atom. We used the conversion relation: 
1 eV/atom = 96.521 kj/mol or equivalently 1 kj/mol = 0.010364 cV/atom. 
Ref. 12. 



s Ref. 54 
h Ref 54 
' Ref 54 
j Ref 52 



m 


Ref. 


[17] 


n 


Ref. 


[17] 


o 


Ref. 


[16] 


P 


Ref. 


[16] 


q 


Ref. 


[15] 


r 


Ref. 


[15] 


s 


Ref. 


[15] 



cubic DOg phase is the stable phase of the synthesized 
Cu3N4>55. 

The odd behavior of the EOS of Ag3N(Fe3N) with 
the existence of two minima (Fig. 1) shows that the 



first minima (to the left) is a metastablc local mini- 
mum that cannot be maintained as the system is de- 
compressed. Ag ions are in the 6g WyckofF posi- 
tions: {x, 0, 0), (0, X, 0), (—X, —X, 0), {~x, 0, (0, ~x, ^) 



6 



Aq^N AqN AgN^ 




O rHOiHO„Li_ DDDDDQcOcOaj DQ'-HfNUUiH^ 
Q<Q— IQ(U-C DQ DQDQ U 

U- CC O 

U 



FIG. 4. (Color online.) Calculated equilibrium properties of the twenty studied phases of silver nitrides. All quantities are 
given relative to the corresponding ones of the fee crystalline elemental silver given in the first row of Table I. The vertical 
dashed lines separate between the difi'erent stoichoimetries. 
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16 18 20 22 24 26 28 30 



P(GPa) 

FIG. 5. (Color online.) Enthalpy H vs. pressure P equation 
of state (EOS) for some AgaN phases in the range where DO9 
— >■ A15, DO9 — >■ DO3 and DO9 LI2 phase transitions occur. 

and (x, x, i); with a; ~ ^ to the left of the potential bar- 
rier (represented by the sharp peak at ~ 18.2 A'^/atom), 
and a; = ^ to the right of the peak. It may be relevant 
to mention here that Cu3N(Fe3N) was found to behave 
in a similar manner^. 

The crossings of the less stable A15, DO3 and LI2 EOS 
curves with the more stable DO9, DO2 and RhFa ones at 
the left side of their equilibrium points reveals pressure- 
induces phase transitions from the latter phases to the 
former. To show this, we plotted the corresponding rela- 
tion between enthalpy H = E{V) +PV and the imposed 
pressure P in Fig. 5. Since DO9, DO2 and RhFs phases 
have identical E{V) curves, the corresponding H{P) 
curves are also identical. Hence, only the H{P) of DO9 is 
displayed in Fig. 5. A point where the enthalpies of two 
phases are equal determine the phase transition pressure 
Pt] and, indeed, the direction of the transition is from 
the higher H to the lower H^^ . As depicted in Fig. 5, 
Pt(D09 -> LI2) = 17.8 GPa, Pt{mg DO3) = 19.5 GPa 
and Pt(D09 ^ A15) = 26.0 GPa. Thus, DO9, DO2 and 
RhFa would not survive behind these Pt's and A15, DO3 
and LI2 are preferred at high pressure. 

Fig. 4 reveals that the AgN group contains the least 
stable phase among all the twenty studied phases: the 
hexagonal B^.. Fig. 2 and Table I show that the sim- 
ple tetragonal structure of cooperite (B17) is the most 
stable phase in this AgN series. In fact, one can see 
from Fig. 4 and Table I that all the considered AgN 
phases possess less binding than their parent Ag(fcc), 
except AgN(B17) which is slightly more stable, with 
0.033 eV/atom lower i?coh- It is interesting to notice 
that AgN(B17) is ^ 0.003 eF/atom more stable than the 
AgsN most stable phases. Moreover, this B17 structure 
was theoretically predicted to be the ground-state struc- 
ture of CuN^, AuN^^ PdN^s and PtN^^. 

Using the full-potential (linearized) augmented plane 
waves plus local orbitals (FP-LAPW+lo) method within 
LDA and within GGA, Kanoun and Said^^ studied the 



E{V) EOS for AgN in the Bl, B2, B3 and B4 struc- 
tures. The equilibrium energies they obtained from the 
E{V) EOS revealed that Bl is the most stable phase, 
and the relative stability they arrived at is in the order 
B1-B3-B4-B2, with a significant difference in total en- 
ergy between B3 and B4 (see Fig. 2(b) in that article). 
Within this subset of structures, the numerical values of 
Ecoh in Table I do have the same order. However, the dif- 
ference between the equilibrium i?coh(B3) and i?coh(B4) 
is only - 0.009 eV , and the E{V) EOS of B3 and B4 
match/overlap over a wide range of volumes around the 
equilibrium point. This discrepancy may be attributed to 
the imphysical/ill-defined measure of stability that Ka- 
noun and Said used, the total energy, while the number 
of the AgN formula units per unit cell in the B4 structure 
differs from that in the others^°. Nevertheless, it may be 
worth mentioning here that AgN(B3) was found to be 
elastically unstable^^. 

In the CuN2 nitrogen-richest phase series, we can see 
from Table I and from Fig. 4 that the phases of this 
group are significantly more stable than all the other 
studied phases, except CI, which is, in contrast, the sec- 
ond least stable among the twenty studied phases, with 
0.017 ey/atom more than AgN(Bk). From Fig. 3, Fig. 
4 and Table I, one can see that in this series, the lower 
the structural symmetry, the more stable is the phase. It 
was found that CuN2 phases have the same trend^. 

Comparing the relative stability of AgsN, AgN and 
AgN2, we find from Table I and from Fig. 4 that AgN2 
in its least symmetric phase, the simple monocline struc- 
ture of CoSb2, is the most energetically stable phase with 
^'coh = —3.699 eF/atom. 

B. Volume per Atom and Lattice Parameters 

The numerical values of the lattice parameters and the 
average equilibrium volume per atom Vb for the twenty 
modifications are presented in Table I. The middle sub- 
window of Fig. 4 depicts the Vb values relative to the 
Ag(fcc). To measure the average of the Ag-Ag bond 
length in the silver nitride, the equilibrium average vol- 
ume per Ag atom (V^"^^), which is simply the ratio of the 
volume the unit cell to the number of Ag atoms in the 
unit cell, is visualized in the same subwindow. 

From the Vb curve in Fig. 4, we can see that, all 
AgN and AgN2 modifications, except the open AgN(Bk) 
phase, decrease Vb; while the AgsN phases tend, on av- 
erage, not to change the number density of the parent 
Ag(Al). 

On the other hand, the ^ curve in Fig. 4 reveals 
that, relative to the elemental Ag and to each other, V^^^ 
tends to increase with the increase in the nitrogen con- 
tent. Thus, in all these nitrides, the introduced N ions 
displace apart the ions of the host lattice causing longer 
Ag-Ag bonds than in the elemental Ag. This cannot be 
seen directly from the Vb values depicted in the same 
figure. 
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For AgN in the Bl, B2, B3 and B4 structures, Ka- 
noun and Said (Ref. 16 described in Sec. Ill A above) 
obtained GGA equilibrium lattice parameters which are 
in very good agreement with ours. Their obtained LDA 
lattice parameter values show the common underestima- 
tion with respect to their and our GGA values (see Table 
I)- 

Gordienko and Zhuravlev^^ studied the structural, me- 
chanical and electronic properties of AgN(Bl) , AgN(B2), 
AgN(B3), AgN2(Cl) and Ag3N(D09) cubic phases. Their 
DFT calculations were based on pseudopotential (PP) 
method within LDA, and on linear combinations of 
atomic orbitals (LCAO) method within both LDA and 
GGA. For comparison, some of their findings are included 
in Table L Within the parameter subspace they consid- 
ered, our GGA values of the a lattice parameter agree 
very well with theirs. On the other hand, although their 
PP a values are closer to the GGA ones (ours and theirs) , 
all their LDA values are less than the GGA ones. This 
confirms the well-known behavior of LDA compared to 
GGA^^"^^. Gordienko and Zhuravlev also found that the 
Ag-Ag interatomic distance increases in the order AgsN- 
AgN-AgN2 . This agrees with the general trend shown in 
Fig. 4, since the V^^^ curve shows an average increase in 
the same direction. 



C. Bulk Modulus and its Pressure Derivative 

Fig. 4 reveals that AgaN phases tend, on average, to 
preserve the Bq value of the parent Ag(Al). Increasing 
the nitrogen content to get AgN phases will increase the 
Bq value of the parent Ag(Al), except in the case of B^. 
While the nitrogen in AgN2 tends to lower the Bq value of 
the parent Ag(Al), the cubic CI phase posses the highest 
Bq value. This could be seen from Fig. 3, where the 
curvature of the £'coh(^) curve of CI is higher compared 
to the shallow minima of the C2, C18 and CoSb2 curves. 

From the definition of the equilibrium bulk modulus 
(Eq. 4), one would expect Bq to increase as i?coii or Vb 
decreases. This is because of the minus sign of the for- 
mer and the inverse proportionality of the latter. That 
is, roughly speaking, the Bq curve should have a mirror 
reflection-like behavior with respect to the i?coh and Vq 
curves. Nevertheless, if E^oh or Vq are increasing and 
the other is decreasing, then the dominant net effect will 
be of the one with the higher change^"*. For example. 
Fig. 4 shows that in going from DO3 to A15, both i?coh 
and Vq increase resulting in a negative change in Bq. In 
going from A15 to DOg, i^coh is decreasing while Vq is in- 
creasing, but, in the end, the latter won the competition 
and lowered the value of Bq. This argument stays true 
throughout the three scries. When there is no significant 
change in both Ecoh and Vq , there is no significant change 
in Bq. This is the case when one goes from CIS to CoSb2. 
A close look at the Bq curve in Fig. 4, reveals that the 
huge decrease in i?coh between CI and C2 defeats the 
relatively small increase in Vq. This is simply because. 



according to Eq. 4, the value of Bq is proportional to the 
absolute change in -Ecoh, while it is far more sensitive to 
any change in Vq because it is proportional to {AVq)~^. 

It is common to measure the pressure dependence of 
Bq by its derivative Bq (Eq. 5). Fig. 4 shows that the 
Bq value of the C2, CIS and CoSb2 phases increases as 
these phases are put under pressure. While the Bq val- 
ues of the rest of the phases shows very low sensitivity to 
pressure and they tend to slightly lower the bulk modu- 
lus, the FcaN phase is the most sensitive phase and tends 
to significantly lower its Bq upon application of pressure. 
This high sensitivity may indicate that the correspond- 
ing minimum on the potential surface is not global, but 
another local minimum as the one at 16.2 A'^/atom (Fig. 

l). 

From the elastic constants they obtained, Gordienko 
and Zhuravlev^^ calculated the corresponding macro- 
scopic bulk moduli (included in Table I). They found 
the highest LDA Bq value for AgN(Bl) among all phases 
they considered, but, in agreement with the present work, 
they obtained the highest GGA Bq value for AgN2(Cl). 
Since LDA relative to GGA overestimates E^oh and thus 
underestimates Vq, each of these two factors (see Subsec- 
tion IIIC) would separately lead to the odd LDA value 
of 219.2 GPa which they obtained. Nevertheless, due 
to this fact, Gordienko and Zhuravlev argued that one 
should consider the LDA and GGA average value of Bq. 



D. Formation Energies 

Formation energies in the present work are used as 
a measure of the relative thermodynamic stabilities of 
the phases under consideration. That is, the lower the 
formation energy, the lower the tendency to dissociate 
back into the constituent components Ag and N2. 

The obtained formation energies Ef oi the twenty re- 
laxed phases are given in Table I and depicted graphically 
in Fig. 4. The latter shows that, relative to each other 
and within each series, the formation energy Ef (defined 
by Eqs. 6 and 7) of the studied phases has the same 
trend as the cohesive energy^^. That is, all phases have 
the same relative stabilities in the Ef space as in the -Ecoh 
space. However, while AgsN phases tend to have equal 
Ecoh as the AgN phases, all AgaN modifications have a 
lower Ef than the AgN ones. Hence, silver nitride is more 
likely to be formed in the former stoichiometry. However, 
all the twenty obtained Ef values are positive; which, in 
principle, means that all these phases are thermodynam- 
ically unstable (endothermic)^^. 

Some of the experimental values of Ef for the syn- 
thesized AgaN phase (which is claimed to be in an 
fee structure) are +199 kJ/mol^^ = 2.062 eF/atom, 
+230 kJ/moli° = 2.3S4 eF/atom, +255 kJ/mol^ = 
2.643 eV/atom and (+314.4+2.5) kJ/mol^ = (3.25853 + 
0.02591) eV; with an average value of 2.587 ± 0.364 eV. 
Among the considered phases in the present work, there 
is only one phase wich has Ef value that fits in this range. 
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the AgN2(Cl). Interestingly, this CI structure has an fee 
underlying Bravia lattice; however, the chemical formula 
differs from that of the synthesized phase. 



E. Electronic Properties 

The DFT(GGA) calculated band diagrams (i.e. ef (k) 
curves) and spin-projected total and orbital resolved (i.e. 
partial) density of states (DOS) of the most stable phases: 
DOg, RhFa, DO2, B17, and C18 are presented in Figs. 6, 
7, 8, 9 and 10, respectively. Spin-projected total density 
of states (TDOS) are shown in sub- figure (b) in each case. 
In all the six considered cases, electrons occupy the spin- 
up and spin-down bands equally, resulting in zero spin- 
polarization density of states: C(e) = n-f{e)~ni{e). Thus, 
it is sufficient only to display spin-up (or spin-down) den- 
sity of states (DOS) and spin-up (or spin-down) band 
diagrams. In order to investigate the details of the elec- 
tronic structure of these phases, energy bands are plotted 
along densely sampled high-symmetry string of neighbor- 
ing k-points. Moreover, to extract information about the 
orbital character of the bands, the Ag(s,p, d) and N(s,p) 
partial DOS are displayed at the same energy scale. 

Fig. 6(a) shows the band structure ef(k) of 
Ag3N(D09). With its valence band maximum (VBM) 
at (i?, —0.086 eV) and its conduction band minimum 
(CBM) at (r, 0.049 eV), Ag3N(D09) presents a semicon- 
ducting character with a narrow indirect band gap Eg 
of 0.134 eV. From sub-figures 6(c) and (d), it is seen 
clearly that the Ag((i)-N(p) mixture in the region from 
— 7.286 eV to —0.086 eV beneath Ep, with two peaks: a 
low density peak around 1.5 eV and a high density peak 
around 4.0 eV steaming mainly from the bands of silver 
d electrons. 

Our obtained PDOS, TDOS and band structure of 
Ag3N(D09) agree qualitatively well with Gordienko and 
Zhuravlev^^; however, using LCAO method within GGA, 
the value of the indirect Eg of Ag3N(D09) they predicted 
is 0.25 eV. 

To the best of our knowledge, there is no experimen- 
tally reported Eg value for Ag3N. However, Tong® pre- 
pared Ag3+a;N samples, and carried out XRD measure- 
ments to confirm the fee symmetry of the prapared sam- 
ples. Using a TB-LMTO code within LDA, Tong then 
calculated the band structure of Ag3N and obtained an 
indirect energy gap of 1.35 eV. Nevertheless, we could 
not figure out the positions of the N ions Tong's model. 

It is a well known drawback of Kohn-Sham DFT-based 
calculations to underestimate the band gap. Thus the 
more demanding GW calculations were carried out, and 
the obtained Eg value will be presented in Sec. Ill F. 

Calculated electronic properties of Ag3N(D02) are dis- 
played in Fig. 8. sub-figure 8(a) shows the energy bands 
ef (k) of Ag3N(D02). With its valence band maximum 
(VBM) at {H, —0.091 eV) and its conduction band mini- 
mum (CBM) at (F, 0.039 eV), Ag3N(D09) presents semi- 
conducting character with a narrow indirect band gap Eg 



of 0.130 eV. From sub-figures 8(c) and (d), one can no- 
tice clearly the Ag(d)-N(p) mixture in the region from 
-7.249 eV to -0.091 eV below Ep, with two peaks: a 
low density peak around —1.3 eV steaming from an al- 
most equal mixture of Ag((i) and N(p), and a high density 
peak around —4.3 eV steming mainly from the bands of 
silver d electrons plus a relatively very low contribution 
from the N(p) states. 

Fig. 7 depicts the band diagram and DOS's of 
Ag3N(RhF3). In contrast to Ag3N(D09) and Ag3N(D02), 
sub-figure 7(a) shows that Ag3N(RhF3) is a semiconduc- 
tor with a narrow direct band gap of 0.129 eV of width 
located at F point. The VBM is at -0.089 eV and the 
CBM is at 0.040 eV. From sub-figures 7(c) and (d), 
one can see the Ag(d)-N(p) mixture is in the region from 
-7.286 eV to -0.089 eV beneath Ep, with two peaks: 
a low density peak around —1.366 eV steaming from an 
almost equal mixture of Ag(c?) and N(p), and a high den- 
sity peak around —4.382 eV steaming mainly from the 
bands of silver d electrons plus a relatively very low con- 
tribution from the N(p) states. 

The relationship between DOg, DO2 and RhF3 struc- 
tures manifests itself in many common features between 
the electronic structure of these three Ag3N nitrides: (i) 
equal Eg of ~ 0.13 eV; (ii) a deep bound band around 
~ —14.6 eV below Ep consists mainly of the N(2s) 
states; (iii) a broad valence band with ~ 7.2 eV of width 
that comes mostly from the 4c? electrons of Ag plus a very 
small contribution from N(2p); and (iv) the relatively low 
TDOS of the conduction bands. 

Energy bands ef (k), total density of states (TDOS) 
and partial (orbital-resolved) density of states (PDOS) 
of AgN(B17) are shown in Figs. 9. It is clear that 
AgN(B17) would be a true metal at its equilibrium. The 
major contribution to the very low TDOS around Fermi 
energy Ep comes from the 2p states of the N atoms as it 
is evident from sub-figure 9(d). Beneath Ep lies a band 
with ~ 7.3 eV of width, in which the main contribution 
is due to the Ag(4(i) states plus a small contribution from 
the N(2p) states. While the N(2s) states dominate the 
deep lowest region around 13.5 eV, the low density un- 
occupied bands stem mainly from the N(2p) states. The 
Fermi surface crosses two partly occupied bands: a lower 
one in the X-M, T-Z-A and T-X-R directions, and a 
higher band in the X-M-T and M-A directions. Thus, 
Ep is not a continuous surface contained entirely within 
the first BZ. 

It may be worth mentioning here that AgN(Bl)^^'^® 
and AgN(B3)^^'^^ phases were also theoretically pre- 
dicted to be metallic. 

Although AgN2(CoSb2) is the most stable phase, but 
the difference in cohesive energy between AgN2(CoSb2) 
and AgN2(C18) is less than 0.02 ey/atom, and we de- 
cided to examine the electronic structure of both phases. 
With Ep crossing the finite TDOS, Fig. 10 shows that 
AgN2(C18) is metallic at K. The orbital resolved 
DOS's reveal that the major contribution to the low 
TDOS at Ep comes from the N(2p) states with tiny 
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FIG. 6. (Color online.) DFT calculated electronic structure for AgsN in the DOg structure: (a) band structure along the 
high-symmetry k-points which are labeled according to Ref. [67]. Their coordinates w.r.t. the reciprocal lattice basis vectors 
are: M(0.5, 0.5, 0.0), r(0.0, 0.0, 0.0), X(0.0, 0.5, 0.0), i?(0.5, 0.5, 0.5); (b) spin-projected total density of states (TDOS); (c) 
partial density of states (PDOS) of Ag(s,p, rf) orbitals in AgsN; and (d) PDOS of N(s,p) orbitals in AgaN,. 
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FIG. 7. (Color online.) DFT calculated electronic structure for AgaN in the RhFa structure: (a) band structure along the 
high-symmetry k-points which are labeled according to Ref. [67]. Their coordinates w.r.t. the reciprocal lattice basis vectors 
are: F(0.5, 0.5, 0.0), 0(0.375,0.625,0.0), 5(0.5,0.75,0.25), Z(0.5, 0.5, 0.5), r(0.0, 0.0, 0.0), L(0.0, 0.5, 0.0), y(0.25, 0.5, -.25), 
E(0.0, 0.5, — .5); (b) spin-projected total density of states (TDOS); (c) partial density of states (PDOS) of Ag(s,p, d) orbitals 
in AgsN; and (d) PDOS of N(s,p) orbitals in AgsN. 




FIG. 8. (Color online.) DFT calculated electronic structure for AgsN in the DO2 structure: (a) band structure along the 
high-symmetry k-points which are labeled according to Ref. [67]. Their coordinates w.r.t. the reciprocal lattice basis vectors 
are: r(0.0, 0.0, 0.0), Af(0.0, 0.0, 0.5), P(0.25, 0.25, 0.25), H(0.5, -.5, 0.5); (b) spin-projected total density of states (TDOS); (c) 
partial density of states (PDOS) of Ag(s, p,ci) orbitals in AgaN; and (d) PDOS of N(s,p) orbitals in AgsN. 
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FIG. 9. (Color online.) DFT calculated electronic structure for AgN in the B17 structure: (a) band structure along the 
high-symmetry k-points which are labeled according to Ref. [67]. Their coordinates w.r.t. the reciprocal lattice basis vectors 
are: X(0.0, 0.5, 0.0), M(0.5, 0.5, 0.0), r(0.0, 0.0, 0.0), ^(0.0, 0.0, 0.5), A(0.5, 0.5, 0.5), 7i(0.0, 0.5, 0.5); (b) spin-projected total 
density of states (TDOS); (c) partial density of states (PDOS) of Ag{s,p,d) orbitals in AgN; and (d) PDOS of N(s,p) orbitals 
in AgN. 




FIG. 10. (Color online.) DFT calculated electronic structure for AgN2 in the C18 structure: (a) band structure along the high- 
symmetry k-points which are labeled according to Ref. [67]. Their coordinates w.r.t. the reciprocal lattice basis vectors are: 
r(0.0, 0.0, 0.0), A' (0.0, 0.5, 0.0), 5'(-.5, 0.5, 0.0), Y{-.5, 0.0, 0.0), ^(0.0, 0.0, 0.5), C/(0.0, 0.5, 0.5), R{-.5, 0.5, 0.5), r(-.5, 0.0, 0.5); 
(b) spin-projected total density of states (TDOS); (c) partial density of states (PDOS) of Ag(s,p, d) orbitals in AgN2; and 
(d) PDOS of N(s,p) orbitals in AgNa. 




FIG. 11. (Color online.) DFT calculated electronic structure for AgNa in the CoSba structure: (a) band structure along 
the high-symmetry k-points which are labeled according to Ref. [67]. Their coordinates w.r.t. the reciprocal lattice basis 
vectors are: r(0.0, 0.0, 0.0), B(-.5, 0.0, 0.0), A(-.5, 0.5, 0.0), £^(-.5, 0.5, 0.5), ^(0.0,0.0,0.5), Y (0.0, 0.5, 0.0), D(-.5, 0.0, 0.5) 
and C(0.0, 0.5, 0.5); (b) spin-projected total density of states (TDOS); (c) partial density of states (PDOS) of Ag(s,p, d) 
orbitals in AgNa; and (d) PDOS of N(s,p) orbitals in AgNa. 
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contributions from the 5s, Ad and 3p states of Ag, re- 
spectively. As one can see from sub-figm-e 10(a), the Ep 
surface crosses the edges of the first Brillouin zone in the 
Z-U-R-S-T-X and T-Z directions. 

The calculated electronic properties of AgN2(CoSb2) 
are displayed in Fig. 11. Band structure, TDOS and 
orbital resolved DOS's have almost the same features as 
the corresponding ones of AgN2(C18). It may be worth 
to mention here that CI phase of AgN2 was also theoret- 
ically predicted to be metalhc^^. 

Compared to the metallic AgN(B17), three new fea- 
tures of these 1:2 nitrides are evident: (i) Deep at 
^ —22.7 eV there is a highly-localized mixture of the 
N(s)-N(p) states. However, the variation in N(2s) en- 
ergy with respect to k is smaller than the variation of 
N(2p) states, resulting in a narrower and higher PDOS. 
(ii) Below the band that is crossed hy Ep there are four 
bands separated by 11.4 eV, ~ 1.6 eV, ~ 0.38 eV 
and 0.28 eV energy gaps, respectively, (iii) The very 
tiny contribution of the N(p) states to the N(2p)-Ag(4(i) 
band. 

A common feature of all the studied cases is that 
Ag(p)-orbitals do not contribute significantly to the hy- 
brid bands. Another common feature is the highly struc- 
tured, intense and narrow series of peaks in the TDOS 
valance band corresponding to the superposition of N(2p) 
and Ag(4(i) states. In their k-space, Ag(4d) energies 
show little variation with respect to k; hence the Van 
Hove singularities-like sharp features. 

To summarize, we have found that the most stable 
phases of AgN and AgN2 are metallic, while those of 
AgsN are semiconductors. A close look at Fig. 9 up to 
Fig. 6 reveals that as the nitrogen to silver ratio increases 
from x = 1 to x = 1/2, the TDOS at Ep decreases; 
and by arriving at a; = 1/3 a gap opens. This finding 
agrees well with Gordienko and Zhuravlev^^. Moreover, 
it may be worth mentioning here that such behavior was 
theoretically predicted to be true for copper nitrides as 
well^'™. 



F. Optical Properties 

Fig. 12 depicts the GW calculated real and imagi- 
nary parts of the frequency-dependent dielectric function 
e-Rphi^) of Ag3N(D09) and the corresponding derived op- 
tical constants. The optical region^^ is shaded in each 
sub- figure. 

The real part £ro(w) (sub-figure 12(a)) shows an up- 
ward trend before 2.3 eV ^ where it reaches its maxi- 
mum value and generally decreases after that. The imag- 
inary part eim(w) (same sub- figure 12(a)) shows an up- 
ward trend before ~ 1.0 eV and it has three main peaks 
located at ~ 2.6 eV in the optical region, ^ 3.3 eV at 
the right edge of the optical region, and at ^ 4.1 eV in 
the UV range. 

Calculated reflectivity R{ijj) and transmitivity T{u)) 
are displayed in sub-figure 12(b). With 0.6 < R{lo) < 0.8, 



it is evident that Ag3N(D0g) is a good reflector, specially 
in the red and the infrared regions. In the visible range, 
the maximum transmitivity T(a;) is at ^ 2.54 eV = 
489 nm, which is at the blue-green edge. 

sub-figure 12(c) depicts the calculated refraction n{uj) 
and extinction K{io) coefficients. As they should, these 
two spectra have, in general, the same qualitative fre- 
quency dependence as the real ei-duj) and the imaginary 
£iin(w) dielectric functions, respectively. 

From the absorption coefficient a (w) spectrum (sub- 
figure 12(d)), it can be seen that Ag3N(D09) starts ab- 
sorbing photons with ^ 0.9 eV energy. Hence, it is clear 
that GWq calculations give a band gap of ~ 0.9 eV, which 
is a significant improvement over the value obtained from 
DFT. The non- vanishing a (w) in the whole optical region 
agrees with the experiment, since it may explain the ob- 
served black color of the synthesized AgaN. 

To the best of our knowledge, the present work is the 
first trial to theoretically investigate the optical proper- 
ties of silver nitride. However, for more accurate optical 
characterization (e.g. more accurate positions and am- 
plitudes of the characteristic peaks), electron- hole exci- 
tations should be calculated. This can be done by eval- 
uating the two-body Green function G2 on the basis of 
our obtained GW one-particle Green function G and QP 
energies, then solving the so-called Bethe-Salpeter equa- 
tion, the equation of motion of 6*2^^. 



IV. CONCLUSIONS 

We have succesfully employed first-principles calcula- 
tion methods to investigate the structural, stability, elec- 
tronic and optical properties of AgaN, AgN and AgN2. 
Within the accuracy of the employed methods, the ob- 
tained structural parameters, EOS, Bq, Bq and electronic 
properties show good agreement with the few avialable 
previous calculations. On the other hand, our obtained 
results show, at least, partial agreement with three ex- 
perimental facts: (i) the lattice parameter of Ag3N(D09) 
is close to the experimentally reported one; (ii) the pos- 
itive formation energies reveals the endothermic (unsta- 
ble) nature of silver nitrides, and (iii) absorption spec- 
trum explains its observed black color. Moreover, the 
present work may be considered as the first trial to the- 
oretically investigate the optical properties of silver ni- 
tride. We hope that some of our obtained results will be 
confirmed in future experimentally and/or theoretically. 
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FIG. 12. (Color online.) The GW calculated frequency-dependent optical spectra of Ag3N(D09): (a) the real £ie(i^) and the 
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